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SUMMARY 


The  results  discussed  in  this  paper  can  be  summarized  as  follows: 

1.  Early  time  (less  than  a  microsecond)  response  of  aerial  lines 

are  apparently  highly  dependent  upon  details  of  the  time- 
dependent  air  conductivity.  Transmission  line  results  do  not 
agree  well  with  finite-difference  calculations  for  this  case. 

2.  For  buried  cables,  early-time  results  from  both  transmission 

line  and  finite-difference  calculations  seem  to  be  in  fair 

agreement.  This  is  probably  due  to  partial  isolation  of  the 
line  from  the  Compton  current  driver  by  the  conducting  ground. 

3.  Late-time  transmission  line  calculations  (for  both  aerial  and 

buried  cables)  indicate  that  large  currents  may  flow  at  late 
times.  This  large  current  is  due  to  the  predicted  slow  decay  of 
the  radial  electric  field  generated  by  a  nuclear  burst,  combined 
with  a  decreasing  line  impedance  as  fields  diffuse  into  the 
ground.  As  a  result,  large  currents  can  persist  for  hundreds  of 
microseconds. 

4.  Quasi-static  analytic  solutions  for  certain  cases  have  also  been 

given.  These  expressions  point  out  the  importance  of  assumed 

loads  at  the  end  of  the  line.  (A  low  impedance  at  each  end  can 

result  in  very  large  late-time  currents.)  Such  solutions  also 

provide  checks  of  the  late-time  accuracy  of  transmission  line 

calculations.  _ . _ _ _ 
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INTRODUCTION 


It  has  been  recognized  for  some  time  that  nuclear  bursts  near  or  at 
the  air-ground  interface  can  create  Compton  current  densities  over  a 
substantial  volume.  This  EMP  source  region  can  result  in  large 
electromagnetic  transients  that  radiate  to  far  distances.  As  a  result, 
complex  EMP  phenomenology  codes  have  been  written  to  calculate  the 
electromagnetic  fields. *“3  Another  concern,  however,  is  the  coupling  of 
such  fields  to  conductors  that  may  exist  inside  the  source  region.  Note 
that  such  conductors  may  include  telephone  or  power  lines  which  will  not 
survive  the  nuclear  blast,  but  are  electrically  connected  to  systems 
outside  the  source  region.  A  large  transient  may  thus  be  transmitted  to  a 
system  of  interest,  even  though  the  conductor  is  destroyed. 

Calculations  of  source  region  coupling  to  long  lines  are  greatly 
complicated  by  the  complex  three-dimensional  geometry,  the  diversity  of 
length  scales  (i.e.,  line  cross-sectional  dimensions  vs.  length),  and  the 
non-linear  air  conductivity  effects.  Many  previous  source  region  coupling 
calculations^-^  have  used  time-domain  transmission  line  models  driven  by 
the  output  of  two-dimensional  EMP  phenomenology  codes.  In  recent  years, 

however,  three-dimensional,  finite-difference  codes®  have  been  developed 
that  can  be  used  to  calculate  at  least  the  early-time  response  of  long 
conductors  inside  the  nuclear  source  region.  Such  techniques  include 

various  non-linear  effects  that  are  missing  in  transmission  line  models. 
This  report  presents  a  comparison  of  such  coupling  calculations. 

EARLY-TIME  FINITE-DIFFERENCE  CALCULATIONS 

Source  region  coupling  is  basically  an  extension  of  surface  burst  EMP 
generation  theory.  Compton  currents  are  used  as  drivers  for  finite- 

difference  approximations  of  the  two  time-dependent  Maxwell  equations, 


while  the  air  ionization  rate  (which  is  proportional  to  the  Compton 
current  at  sea  level  air  densities)  is  used  in  simplified  air  chemistry 
reaction-rate  equations  to  calculate  the  time-dependent  electron  and  ion 
densities.  These  densities  are  then  combined  with  field-dependent 
expressions  for  the  mobility  to  obtain  a  time-varying  conductivity,  which 
is  used  in  Maxwell's  equations.  Details  of  this  theory  have  been 
discussed  elsewhere, 1  and  will  thus  not  be  repeated  here.  Note, 
however,  that  such  calculations  are  non-linear  primarily  due  to  the  field- 
dependent  mobility. 

Unfortunately,  adding  even  a  radial-directed  wire  to  a  surface  burst 
EMP  calculation  destroys  the  usually  assumed  2-D  symmetry  and  greatly 
increases  the  numerical  complexity.  At  early  times  (less  than  about  a 
microsecond),  however,  the  high  air  conductivity  makes  the  problem 
primarily  a  local  one  in  that  the  Compton  currents  in  a  small  volume 
around  the  conductors  generate  the  fields  which  are  responsible  for 
currents  on  the  conductors.  In  this  time  regime,  3-D  finite-difference 
time-domain  (FDTD)  calculations  are  possible,  although  special  techniques 
are  required  to  reduce  computer  costs  to  reasonable  values. 

As  an  example,  the  PRES-3D  computer  code®  was  used  to  calculate 
SREMP  coupling  to  the  aerial  wire  geometry  shown  in  Figure  1.  The  weapon 
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Figure  1.  Geometry  for  SREMP  coupling  to 
a  wire  over  a  ground  plane. 


is  assumed  to  be  detonated  approximately  350  m  away  from  the  cable  and  on 
the  axis  of  the  cable.  The  cable  length  was  selected  to  be  long  enough  to 
allow  observation  of  current  propagation,  but  was  limited  by  computer  cost 
considerations.  The  cable  was  placed  10  meters  above  a  perfectly  con¬ 
ducting  ground.  The  gamma-flux  was  assumed  to  be  parallel  to  the  cable. 
At  a  point  z  meters  from  the  weapon  and  t  seconds  after  detonation,  the 
Compton  current  Jz(z,t)  and  air  ionization  rate  S(z,t)  in  the  vicinity 
of  the  cable  were  prescribed  as 


Jz(z,t)  =  3.4  x  io10  £XP(.Z?.Z120).  f(t-2/c)  (A/m2) 


$ 


S(z,t)  =  8.2  x  i()22  j(Zjt)  (ion  pairs/m^/s)^  .  (2) 

The  time  history,  f(t),  and  the  peak  Compton  current  versus  axial  position 
are  displayed  in  Figure  2.  The  rise  time  of  the  current  density  and  the 
cable  dimensions  force  grid  sizes  for  the  finite-difference  solution  to  be 
on  the  order  of  10  cm  to  1  m.  For  a  cable  about  1  km  in  length,  many 
millions  of  such  points  would  be  required  to  solve  the  finite-difference 
equations  in  the  entire  volume.  This  is  not  practical  on  existing 
computers.  If  the  cylindrical  wire  is  perfectly  conducting,  any  induced 
fields  will  be  propagated  along  the  wire  axis  with  a  velocity  determined 
solely  by  the  external  medium.  Thus,  to  examine  the  prompt  fields  on  the 
wire,  it  is  not  necessary  to  model  the  entire  cable  for  all  times,  but 
only  a  small  fraction  of  the  cable  length  which  encompasses  the  Compton 
current  pulse  and  which  moves  along  the  wire  with  the  pulse. 


i  (*) 


Figure  2.  The  noraalized  Coaptou  current  tiae  history  f(t)  and 
the  peak  Coapton  current  density  vs.  axial  position. 


One  thus  has  a  "moving  window"  simulation  volume  which  propagates 
down  the  wire  at  the  velocity  of  the  driver.  Such  a  technique  enables  one 
to  Include  many  of  the  advantages  of  a  retarded  time  calculation  in  a  real 
time  computer  code.  Also,  using  this  method,  and  one  plane  of  symmetry 
through  the  wire,  reduces  the  number  of  grid  points  to  a  manageable  level. 

For  this  problem,  the  moving  simulation  volume  was  75  m  along  the 
axis  of  the  wire.  This  problem  has  a  plane  of  symmetry  in  a  vertical 


plane  through  the  wire  center,  which  reduced  the  number  of  required  mesh 
points  by  a  factor  of  2.  The  zoning  along  the  y  axis  was  carried  out  to 
31  m  from  wire  center.  Along  x  (vertical),  there  was  10  m  between  wire 
and  ground  and  21  m  between  the  wire  and  the  top  of  the  simulation 
volume.  All  boundaries  of  the  simulation  volume  had  perfectly  conducting 
boundary  conditions.  The  use  of  this  boundary  condition  for  the  ground 
was  intentional  to  allow  compa  'sons  between  this  calculation  and  other 
models.  It  is  expected  that  the  highly  conducting  air  during  the  gamma 
pulse  will  not  allow  the  other  artificial  boundary  conditions  to  affect 
the  results  for  our  short  timeframe  of  interest  in  this  calculation.  The 
moving  simulation  volume  used  67,000  zones.  The  zoning  along  the  z  axis 
necessary  to  resolve  the  propagating  source  was  a  constant  1.5  m.  Along 
the  x  and  y  axes,  the  zoning  was  11  cm  in  the  vicinity  of  the  wire  and  was 
gradually  extended  to  several  meters  at  the  extreme  of  the  simulation 
volume. 

The  Compton  currents  were  defined  at  the  mesh  points  usir.g  the 
analytical  formula  of  Equation  (1).  The  ionization  rate  of  Equation  (2) 
wds  used  as  the  source  term  for  the  three-species  air-chemistry 
equations.  The  solutions  between  the  field  equations  and  the  air- 
chemistry  equations  were  done  self-consistently — i.e.,  the  air-chemistry 
parameters  and  conductivity  are  field-dependent  and,  similarly,  the  field 
solutions  depend  strongly  on  the  conductivity. 

To  anticipate  some  of  the  results  of  the  simulation,  it  was  observed 
that  neither  the  perfectly  conducting  boundary  conditions  nor  the  large 
cells  at  the  end  points  appeared  to  have  any  influence  on  the  current 
carried  by  the  wire.  This  is  no  doubt  due  to  the  very  long  path  that  any 
EM  wave  must  transverse  through  an  ionized  conducting  medium  from  a 
boundary  to  the  wire. 


Figure  3  displays  the  peak  value  of  the  cable  current  as  a  function 
of  axial  position  along  the  wire.  Time  histories  of  the  current  at 
several  points  along  the  wire  are  also  displayed  in  this  Figure.  The 
current  waveform  rises  rapidly  (~40  nsec),  then  tends  to  be  relatively 
constant  for  about  100  nsec,  and  then  starts  to  decay. 
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Figure  3.  Peak  value  of  the  cable  current  vs.  axial  position  and 
tiae  variation  of  the  current  on  the  cable  at  a  few  axial 
positions  (arrows  on  tiae  histories  indicate  window  tiaes 
at  that  z  location). 


Figure  4  displays  the  peak  value  of  the  air  conductivity  versus  axial 
position  for  a  point  near  the  cable  and  a  point  between  the  cable  and 

ground.  Note  that  the  conductivity  is  significantly  higher  near  the  cable 

because,  as  we  approach  the  cable  (or  the  ground),  the  boundary  conditions 
require  that  the  tangential  electric  field  be  zero  on  the  cable  and 

ground.  The  lower  field  values  near  the  ground  and  the  wire  give  much 


higher  electron  mobilities  and,  hence,  a  higher  conductivity.  (Zero  water 
vapor  content  was  assumed  in  these  calculations.) 


•  !•] 


Figure  4.  Conductivity  vs.  axial  position  near  the 
cable  and  between  cable  and  ground. 

Note  that  the  aerial  cable  of  this  calculation  is  10  centimeters  in 
diameter  ard  is  uninsulated.  The  high  conductivity  near  its  surface  thus 
makes  it  possible  for  large  currents  to  flow  on  or  off  the  cable.  The 
situation  might  be  even  more  complicated  if  the  wire  were  very  thin 
because  large  perpendicular  electric  fields  would  exist  at  the  surface  of 
the  wire,  perhaps  reducing  the  conductivity  there.  At  a  few  centimeters 
from  the  surface,  however,  fields  would  be  small,  giving  an  enhanced 
conductivity.  In  any  case,  the  air  conductivity  in  any  plane  perpen¬ 
dicular  to  the  axis  of  the  long  conductor  is  apparently  not  uniform  (a 
fact  often  assumed  in  simpler  models). 


As  a  second  example,  finite-difference  time-domain  techniques  were 
used  to  calculate  the  response  of  an  insulated,  buried  cable.  The  cable 
was  10  meters  in  length  and  was  buried  one  meter  below  the  ground-air 
interface.  The  ground  conductivity  was  assumed  to  be  10“^  mho/m  and  the 
dielectric  constant  was  6e0.  The  conducting  cable  was  1/4  inch  in 
diameter,  and,  for  these  calculations  the  source  term  was  a  prescribed 
electric  field  along  the  air-ground  interface  above  the  cable.  The  time 
history  of  this  field  was  triangular,  rising  to  a  peak  of  100  kV/m  in  5 
nsec,  and  then  returning  to  zero  at  300  nsec.  The  starting  time  of  this 
prescribed  surface  field  was  phased  at  the  speed  of  light  along  the 
cable.  Some  sample  results  for  this  example  are  shown  in  Figure  5.  As 
before,  calculations  are  limited  to  early  times  by  the  expenses  of 
computer  time. 

TRANSMISSION  LINE  CALCULATIONS 

Transmission  line  raodeling^-13  a  fairly  well  established  method 

for  calculating  EMP  coupling  to  long  conductors  that  are  outside  of  the 
nuclear  source  region.  Such  methods  are  also  used  for  coupling  inside  the 
source  region,  but  the  adequacy  of  such  methods  is  rather  uncertain  under 
such  conditions  (especially  at  early  times). 

Transmission  line  techniques  have  the  great  advantage  of  being 
capable  of  calculating  late-time  response  of  very  long  conductors  at  a 
reasonable  cost.  Such  late-time  response  calculations  are  believed  to  be 
very  important  because  EMP  generation  models  indicate  that  substantial 
late-time  electric  fields  are  produced  by  a  surface-level  nuclear  burst. 
The  apparent  result  is  very  large  amounts  of  late-time  (~1  millisecond)  EM 
energy  coupled  to  long  lines. 


Figure  5.  Currents  on  a  burled  cable  at  early  tines  as  calculated  by  FDTD  code. 


The  time-dependent  transmission  line  equations  can  be  written*®  as 


l-5F+  ri 


v  -  ** 

vd  17 


(3) 


and 


3V 

C17  +  GV 


_3I 

3z 


(4) 


where  I  and  V  are  the  current  and  voltage  along  the  line;  L,  C,  R  and  G 
are  respectively  the  inductance,  capacitance,  resistance,  and  conductance 
per-unit-length  of  the  line  and  Vq  is  the  per-unit-length  voltage  source 
driver.  The  driver,  VD,  is  usually  assumed  equal  to  the  electrical 
field  along  the  wire  if  the  wire  was  absent  (i.e.,  the  incident  electric 
field). 


Note  that  in  many  special  cases,  one  can  obtain  analytic  solutions 
for  the  current  and  voltages  along  the  line.  Numerical  solutions  are  also 
straightforward  and  with  implicit-differencing  techniques,  solutions  can  be 
obtained  over  long  time  periods.  The  primary  problem  in  using  trans¬ 
mission  line  theory  is  thus  determining  the  proper  expressions  for  four 
line  parameters  L,  C,  R  and  G. 


In  some  cases,  appropriate  line  parameters  are  well  known.**  If 
one  assumes  an  aerial  conductor  of  diameter  d  at  a  height  h  above  a 
perfectly  conducting  ground  plane,  then 

.  uo  -1  /  2h\ 

L  -  ^  ^sh  ( j-J 


(3) 


while  if  the  region  between  the  wire  and  the  ground  has  a  uniform 
conductivity,  a,  then 


(7) 


For  this  example,  currents  will  flow  in  one  direction  along  the 
aerial  wire  and  return  through  the  ground  plane,  with  some  leakage  through 
the  intervening  conducting  media.  For  a  perfectly  conducting  ground,  the 
resistance  per  unit  length  is  zero. 

One  approach  is  to  simply  use  these  expressions  with  a  time-dependent 
conductivity  in  Equation  (7),  thus  giving  a  time-dependent  conductance. 
The  conductivity  to  be  used  might  be  that  calculated  by  a  2-D  EMP  code 
which  does  not  include  the  wire.  (Such  a  calculation  also  provides  the 
per  unit  length  voltage  source,  Vq.) 

Such  an  approach  does  not  include  the  fact  that  at  early  times  in  the 
source  region,  air  conductivities  become  large  enough  to  isolate  an  aerial 
conductor  from  the  local  ground.  Return  currents  will  then  flow  through 
the  conducting  air,  rather  than  on  the  ground  plane. 

An  empirical  method  for  attempting  to  include  this  effect^  is  to 
define  a  time-dependent  skin  depth  in  the  air,  5(t),  where 

5<t>  ■  Vw  •  <8) 

When  <j(t)<h,  it  can  be  assumed  that  return  currents  flow  through  the  air 
within  a  distance  <5  of  the  wire. 


One  then  uses  the  coaxial  formulas ^ 


WtV 


2iro(t) 

mn 


o(t)ir62(t) 


On  the  other  hand,  when  the  skin  depth  becomes  larger  than  the  wire 
height  above  the  ground,  the  previous  expressions  for  the  line  parameters 
are  applicable. 

Some  representative  results  of  a  transmission  line  calculation  using 
these  expressions  are  shown  in  Figure  6.  In  this  case,  the  assumed 
geometry  is  identical  to  that  previously  used  in  the  finite-difference 
calculations  (see  Figure  1)  while  the  time-dependent  source  term,  Vp, 
and  conductivity,  o(t)  used  in  the  transmission  line  calculation  were 
taken  from  the  finite-difference  results  (at  a  point  distant  from  the 
aerial  conductor). 


V?! 


Time  (us) 


Figure  6.  Early-time  current  wave fo rue  as  a  function  of  position 
on  an  aerial  cable  as  predicted  by  a  transmission  line 
calculation.  Geometry  is  the  same  as  in  Figure  1. 

A  comparison  of  these  transmission  line  results  with  those  previously 
obtained  using  finite-difference  techniques  shows  little  agreement. 
Transmission  line  currents  rise  relatively  slowly  but  continually  over  the 
time  period  calculated,  compared  to  the  early  time  peaks  seen  in  the 
finite-difference  calculations.  However,  one  would  somewhat  expect  the 
transmission  line  calculations  to  be  incorrect  at  early  times  because  many 
approximations  and  assumptions  inherent  to  such  models  are  clearly 
invalid.  For  example,  the  air  conductivity  in  any  plane  perpendicular  to 
the  conductor  was  assumed  to  be  independent  of  location,  an  assumption 
that  is  demonstrated  to  be  incorrect  by  the  finite-difference  results  of 
Figure  4. 


Transmission  line  techniques  may  be  more  applicable  for  a  buried 
conductor,  however,  where  the  local  conductivity  is  fairly  uniform  and 
does  not  vary  with  time.  Fields  generated  by  Compton  currents  in  the  air 
diffuse  into  the  ground  to  drive  the  cable,  while  return  currents  must 
also  flow  through  the  conducting  soil.  Again  one  can  use  the  concept  of  a 
time-dependent  skin-depth 


V*) 


where  <Jg  is  the  ground  conductivity.  As  before,  this  results  in  time- 
dependent  transmission  line  parameters  (see  Equations  (9)  through  (12)). 
However,  L,  C  and  G  have  only  a  logarithmic  time  dependence,  while  the 
resistance  per  unit  length,  R,  is  approximated  by 


"Og  <S|(t)  2irt 


The  resistance  of  the  current  return  path  through  the  ground  thus 
decreases  linearly  with  time. 

Results  of  a  transmission  line  model  of  coupling  to  a  10  meter  long 
cable  buried  one  meter  deep  are  shown  in  Figure  7.  Problem  parameters  are 
the  same  as  previously  discussed  with  the  comparable  finite-difference 
calculation.  A  comparison  of  these  results  with  those  previously  shown  in 
Figure  5  indicates  a  reasonable  agreement  between  the  simple  transmission 
line  calculation  and  finite-difference  results.  Predicted  waveforms  have 
quite  similar  shapes,  although  magnitudes  differ  by  about  a  factor  of  two. 


Consider  now  a  long  buried  conductor  and  its  SREMP  response  at  late 
times  (i.e.,  times  of  roughly  1  millisecond).  Finite-difference  time¬ 
stepping  calculations  cannot  practically  be  carried  out  to  such  late 
times,  but  transmission  line  methods  can  be.  Fortunately,  such  results  at 
late  times  can  be  checked  against  certain  analytic  calculations  or  circuit 
models. 

As  an  example,  consider  a  conductor  (either  buried  or  in  the  air)  of 
length  i  (where  i  may  be  of  the  order  of  1-2  kilometers).  If  late-time 
SREMP  fields  are  slowly  changing,  skin  depths  and  wavelengths  are  greater 
than  or  equal  to  cable  lengths.  Equation  (14)  indicates  that  the 
resistance  per  unit  length  of  the  line  is  small,  as  is  the  inductive 
impedance  at  late  times.  For  an  insulated  line,  the  late-time  current 
then  depends  primarily  on  the  loads  at  the  ends  of  the  conductor.  If  E0 
is  the  average  incident  electric  field  along  the  conductor,  and  Rj  and 
R2  are  resistive  loads  at  each  end,  then  the  late-time  current  is 
approximated  by 


’  RX  +  R2  +  R* 

as  is  readily  seen  by  a  simple  equivalent  circuit  model.  If  such  a  line 
runs  radially  outward  from  a  nuclear  fireball,  the  load  Rj  may  be  quite 
small  (i.e.,  the  line  is  shorted  to  ground  by  the  conducting  fireball). 
The  current  at  late  times  is  thus  primarily  limited  by  the  load  at  the  far 
end,  R2.  EMP  environment  codes  indicate  that  substantial  electric 
fields  may  exist  at  late  times.  For  example,  an  E0  of  even  100  v/m  over 
a  kilometer  of  cable  gives  a  voltage  source  of  10^  volts  and  if  R2  is 
about  100  ohms,  late-time  currents  of  1000  amps  will  be  generated. 


If  a  cable  is  uninsulated,  late-time  results  are  slightly  more 
complicated  because  current  can  leak  off  the  cable  along  its  entire 
length.  However,  in  the  static  limit,  when  a  uniform  driver  is  applied  to 
an  uninsulated  transmission  line,  one  can  show  that  the  spatial  current  is 
given  by 


I(z)  -  |a  exp( /RGz)  -  B  exp(-/RGz) j 


(16) 


where  A  and  B  are  constants  which  depend  upon  the  loads  at  each  end  of  the 
line  (i.e.,  A  and  B  are  determined  by  the  boundary  conditions).  For 
example,  with  resistive  loads  Rj  and  R2,  where  Ri=0,  then  A*  -B  and 
one  can  show  that 


A 


EqR2/R 
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(17) 


where  d  -  *TlG£  and  l  is  the  length  of  the  line.  This  solution  indicates 
that  there  is  a  characteristic  length  at  late  times  given  by  (RG)-*/^ 

Late-time  transmission  line  calculations  can  thus  be  checked  against 
such  specific  analytic  solutions.  As  an  example,  consider  a  2000  meter 
long  cable,  buried  in  10“^  mho/m  soil,  driven  by  the  field 


Eq(t)  -  100( t/T)“1/2  v/m 


(18) 


where 


t  ■  retarded  time  ■  t  -  z/c 


Transmission  line  calculations  of  the  resulting  cable  currents  for 
insulated  and  uninsulated  lines  are  shown  in  Figures  8  and  9.  A 

comparison  of  these  results  at  late-times  with  the  analytic  expressions  in 
Equations  (15)  and  (16)  indicate  that  the  transmission  line  results  do 
indeed  approach  the  analytic  expressions  for  the  static  limits. 

It  is  also  interesting  to  note  that  currents  along  most  of  the  length 
of  the  uninsulated  wire  exceed  those  on  the  insulated  wire.  At  the  far 
end  (at  R2),  however,  the  current  is  smaller  for  the  uninsulated  case 
than  for  the  insulated  line. 


(  I  .  .  .  '  ♦  •  '  •  1  • 

TIME  (ms) 


Line  Parameters:  Driver: 

l  -  2000  m  EQ( T)  =  100  (t/T)-1/2  V/m 

=  1  ohm 

R,  *  10  ohms  T  =  10-2  s 

*- 

Solid  line  -  transmission  line  results 

*  -  prediction  of  quasi-static  analytic  theory 

Figure  8.  Current  at  the  center  of  a  2000  meter  long  insulated 
cable  (current  at  other  locations  along  cable  is 
virtually  identical). 


Figure  9.  Currents  on  a  2000-aeter  buried  cable  in  good  electrical 
contact  with  the  ground  driven  by  a  distributed  source. 
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